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CN Abstract 

q We examine four dimensional magnetically charged extremal black holes in certain non-linear U(l) gauge 

theories coupled to two derivative gravity. For a given coupling, one can tune the magnetic charge (or vice 

versa) so that the curvature singularity at the centre of the space-time is cancelled. Since these solutions 

have a horizon but no singularity, they have been called regular black holes. Contrary to recent claims in 

the literature, we find that the entropy function formalism reproduces the near horizon geometry and gives 

(— i the correct entropy for these objects. 
+-> 

1 Introduction 

i-H The Penrose cosmic censorship hypothesis states that, if singularities predicted by General Relativity occur 
^ in nature, they must be dressed by event horizons [1]. Behind the veil of an event horizon, there is no causal 
j-^ contact from the interior to the exterior of a black hole, so the pathologies occurring at the singular region 
C*") can have no influence on an external observer. However, the converse of the hypothesis is apparently not 
true — a horizon does not necessarily hide a singularity. Solutions with a horizon but no singularity have 
,—l been called regular black holes. 

^ The holographic principle, [2, 3], states that the number of degrees of freedom describing the black hole 
<0 is bounded by the area of the horizon. A stronger statement is that degrees of freedom living on the horizon 
;> can describe the physics of the interior completely. While the holographic principle is essentially a proposed 
feature of quantum gravity, one might wonder whether having a classically regular or singular solution has 
any quantitative or qualitative effect on the entropy and the physics at the horizon. 

Earlier work on regular black hole models can be found in [4-8] These regular solutions are referred to as 
Bardeen black holes [9] . In addition, regular black hole solutions to Einstein equations with various physical 
sources were reported in [10] and [11]. Among known regular black hole solutions, are the solutions to the 
coupled equations of nonlinear electrodynamics and general relativity found by Ayon-Beato and Garcia [12] 
and by Bronnikov [13]. The latter describes magnetically charged black hole, and provides an interesting 
example of the system that could be both regular and extremal. In this note we are specially interested in 
the near horizon geometry of an extremal magnetically charged black hole non- linearly coupled to a U{1) 
gauge field. For a given magnetic charge, one can tune the non-linear coupling so that the solution is regular. 

In [14], Matyjasek found the near horizon, AdSz x S 2 geometry of a particular magnetically charged 
extremal black hole. The entropy function formalism of Sen [15, 16], is particularity useful for discussing 
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the entropy of extremal black holes especially when non-linear or high derivative terms make a full analysis 
difficult. Since, the formalism is equivalent to solving Einsteins equations for the near horizon region, a priori, 
and assuming the near horizon geometry decouples, one would expect to be able to reproduce the results of 
[14] using Sen's approach. This issue has been studied recently, [17, 18], and authors reported that, even at 
the level of two derivative gravity, the entropy function approach does not lead to the correct Bekenstein- 
Hawking. To account for this discrepancy, they claim that the entropy function approach is sensitive to 
whether the nature of the central region of the black hole is regular (linear) or singular (nonlinear). 

Contrary to the claims of, [17, 18], in this note we find that a straight forward application of the entropy 
function formalism reproduces the results of [14]. The equation of motion derived from extremizing the 
entropy function are exactly the same as equation of motion at horizon found by extremizing the action, 
since the entropy function (up to Legendre transformation) is the Lagrangian at the horizon. The fact that 
the entropy is the value of entropy function at its extremum is derived from the Wald entropy formula, 
using the near horizon symmetries. Both of these results, just coming from careful consideration of the near 
horizon symmetries and have nothing to do with the regularity of the solution inside the horizon. Further 
more, we find that by varying the non-linear coupling, the regular solution can be smoothly connected to 
the extremal Reisner-Nordstrom solution of Einstein-Maxwell theory. 

The paper is organised as follows. In section 2 we review a particular regular black hole solution of 
interest. Then, in section 3, we review the entropy function formalism and apply it to Einstein gravity 
coupled to non-linear electrodynamics. In section 4 we consider the special case of the formalism applied 
to a regular black hole solution. Finally we end with the conclusion in section 5, having relegated some 
technical details about the large charge, small coupling expansion of the entropy to appendix A. 



2 Regular Black holes 

In this section we review a magnetically charged regular black hole solution of Einstein gravity coupled to 
non-linear electrodynamics and its extremal limit [12-14, 19], mainly following [14] with slightly different 
notation. 

We consider an action given by, 

S= lLJ d4x ^ R ~ £ H^ 2 )) , (2-1) 
where F 2 = F^F^ and the non-linear U(l) gauge field Lagrangian, Cf, is, 

£ F = F 2 cosh- 2 ((A 2 F 2 /2) 1/4N ) • (2.2) 



The equations of motion corresponding to the metric and gauge field and the Bianchi identity are, 

Rpu - 2 g ^R = Q(p2^ F L lxF v X ~ ^Fdiiv , (2.3) 

d [tl F aP] = . (2.5) 

For a magnetically charged black hole, the equation of motion for the gauge field and the Bianchi identity 
can be solved by, 

F e4> = Psme, (2.6) 
where P is the magnetic charge of the black hole. A static, spherically symmetric ansatz for the metric: 



ds 2 = -a 2 {r)dt 2 + + r 2 d£l\ , (2.7) 

can solve Einstein equations with, 



dr 2 
a 2 (r) 



9/ s 2m(r) 

a 2 r = 1 — , 2i 

r 

lr The coupling a — \/A is commonly used in the literature. With out loss of generality, we can take A > 0. 
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where, 

\P\ lA/Pl 1 / 2 

m M = m 0O 1 1 tanh 1 1 , . (2.9) 

V ' 2|A/P| 1 /2 r /\P\ V ; 

The parameter, moo, is an integration constant which can be fixed by employing the boundary condition 
m(oo) = M, where M is the black hole mass. Moreover demanding of the regularity of the line element as 
r — ► 0, yields, 

M = wh>' (2 - 10) 

and consequently, m(r) reads, 

m(r) = M (l — tanh = — ^j— ( 1 - tanh ] . (2.11) 



2Mr J 2|A/P| 1 /2 I r /\P\ 

The location of the inner and outer horizons, r±, which are given by equation a{r) = 0, can be expressed in 
terms of the real branches of the Lambert function, Wi(x), as follows, 

2 2 

r + p r_ p 

M = ~Wo(-eP 2 /V/ 4 ) _ p 2/ 4 ' M = ~W_i(-eP 2 /V/4) -p 2 /4 ' 

where, p = P/M, is the magnetic charge-to-mass ratio. The Lambert function 2 , is defined by the formula, 

e w{x) W(x) = x . (2.13) 

This function has two real branches, called Wo and W-\, with the branch point at x = —1/e. Since the 
value of the principal branch of the Lambert function, Wq, at 1/e, plays an important role in our discussion, 
we define wo = Wo(l/e). 

1/2 

When p = p ex t = 2w , r+ = r_, and the two horizons merge into a degenerate horizon giving an 
extremal solution. Since we will be considering the near horizon geometry, we will eliminate the mass from 
our formulae as it is defined asymptotically. Using (2.10) we can express the condition for extremality and 
regularity, p ext = 2u;q /2 , as, 

p = w . (2.14) 

In other words for an extremal black hole to be regular we must tune the charge to coupling ratio to a 
particular value. A generic extremal, but not necessarily regular, solution to (2.1) will still have a degenerate 
horizon but presumably with a different charge to coupling ratio. 
One can write the near horizon limit, found by [14] as: 



ds 2 = v\ 



with, 



-p 2 dt 2 + + v 2 {dO 2 + sin 2 6># 2 ) . (2.15) 



v 2 = j^\ 2 P 2 - 0.68P 2 , (2.16) 

(1 + Wq) z 



vi = jv^sP 2 - 1-07P 2 , (2.17) 
(1 + w ) 6 

V 2 1 



* ! 1 

and the Bekenstein-Hawking entropy is, 



1 + w ) « 0.64 , (2.18) 



(1 + wo) 2 

2 See [20] for a nice review of the properties of the Lambert function 



S BH = \A = 7rv 2 = r^^P 2 . (2.19) 
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3 Entropy function Analysis 



In this section we briefly review the entropy function formalism of Sen [15, 16] and subsequently apply it to 
magnetically charged extremal solutions of (2.1). 

Assuming a gauge and diffeomorphism invariant Lagrangian and a near horizon AdS2 x S 2 geometry, the 
entropy function is defined as the Legendre transform, with respect to the electric charges, of the reduced 
Lagrangian evaluated at the horizon: 



£(u, v, Q, P) = 2vr e'Qi - f(u, v, e, P) = 2vr e l Q t - / d0dipV^GC , (3.1) 




where e l are the electric fields, Qi = df/de l , are the electric charges conjugate to the electric field, u are the 
values of the scalar moduli at the horizon, and v\, v 2 are the sizes of the AdS2 and S 2 , respectively. The 
near horizon equations of motion for a black hole carrying electric charges Q and magnetic charges P, are 
equivalent to the extremisation of £ with respect to u, v and e: 

d£ d£ d£ , . 

^- = 0, 7^ = 0, ^ = 0. 3.2 
ou ovi oe 

Furthermore, the Wald entropy associated with the black hole is given by £ at the extremum (3.2). If 
£ has no flat directions, then the extremization of £ determines u, Vi, and e, in terms of Q and P. The 
extremal value of the Wald entropy, S = £(Q, P)\ e xtr, is independent of the asymptotic values of the scalar 
fields. This neatly demonstrates the attractor mechanism, [21-23], with out requiring supersymmetry [24]. 
The formalism can even be extended to rotating black holes which have less near horizon symmetry [25]. 
However, since it only involves the near horizon geometry, a weakness of the formalism is that one implicitly 
assumes that the full solution exists, which is not always the case [26]. 

We now specialise our discussion to the case of interest. Since the regular black hole solution has an 
extremal limit, one can use the entropy function formalism to find the near horizon geometry and the 
entropy. We take the near horizon AdS2 x S 2 metric to be give by (2.15). From the definition (3.1), using 
the Lagrangian (2.1), the entropy function evaluates to, 

£ = vr (v 2 - vi + \viv 2 C F (2P 2 /v 2 2 )) . (3.3) 

By extremizing this entropy function with respect to v\ and V2, we find following equations, 

= -1 + l 2 v 2 C F {2P 2 /v 2 ) , (3.4) 

= l + \v 1 ^[v2£ F (2P 2 /v 2 )] . (3.5) 

Substituting (3.4) into (3.3) gives, 

£ = nv 2 = \A , (3.6) 

which is just the Bekenstein-Hawking entropy. This result, which is independent of the form of Cp, is to 
be expected, since, in the absence of higher derivative terms, the Bekenstein-Hawking and Wald entropies 
coincide. 

Now, the equations of motion allow us to determine vi and V2 in terms of P and the coupling A. The 
first equation, (3.4), determines V2, and consequently the entropy, in terms of P (and A). Having found V2, 
(3.5) allows us to determine v\ in terms of V2- Consequently, we see that extremising the entropy function 
completely determines the entropy and near horizon geometry in terms of P (and A). 

We now consider explicitly finding v± and t>2 for a particular Lagrangian. Using the Lagrangian (2.2), 
(3.4) and (3.5) become, 

= -1 + (P 2 /v 2 ) cosh- 2 (VAPM) , (3.7) 
= 1 - Vl {P/v 2 f cosh~ 2 {^\P/v 2 ) 

+ V! V\ {P/V 2 ) 5/2 COsh^ 3 ( y/ \P / V 2 ) Smh ( yj XP / V 2 ) , (3.8) 
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which agrees with the near horizon equations of motion found directly in [14]. 
To solve (3.7), it is convenient to rewrite it as, 

cosh£ = 7£, (3.9) 

where, 

£ = ^APM, 7 = (A/P)- 1/2 . (3.10) 

One can then graphically solve (3.7) by finding the intersection of cosh£ and 7^ for various values of 7. We 
illustrate this procedure in figure 1. It is not hard to see that as we increase 7, there are either zero, one 
or two solutions to (3.9). One can also see from figure 1, that, as X/P — * (i.e. 7 — * 00) , the two possible 
values for £ are, 

ei^o-jo- (3.11) 

Notice that, since cosh 2; > 1, (3.7) also implies, 

v 2 < P 2 . (3.12) 



cosh(£) 




Figure 1: This figure illustrates the graphical solution of (3.9) which is given by the intersection of cosh£ and 7^. 
As we increase the gradient, 7 = y/P/X, one obtains either no solutions, a tangential point (denoted by a red square 
above) or two solutions. The first point of a double intersection, labelled with a brown dot, corresponds to a point on 
what we call, for reasons that will be clear later, the large branch and the second intersection, labelled with a green 
triangle, is on the small branch. 



Now, we can define a (multi- valued) function, ^(x), by, 



cosh.F(x) 

so that we can formally write down a solution to (3.9) as, 
Then letting, 

Q(x) 

and using and (3.10), we can write, 



X (3.13) 



1 



cosh {T{x)) 



v 2 



(3.13) 



(3.14) 



(3.15) 



(3.16) 
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which is of the generic form expected by dimensional analysis. Since (3.9) may have two solutions, T and 
Q both have two branches. Substituting (3.11) into (3.15) we find that, in the limit that the non-linear 
coupling goes to zero (or the charge becomes very large) , 



0(0) = l/cosh 2 (^(0)) 



(small branch, Qg) 

1 (large branch, Ql) 



(3.17) 



We call the two branches of Q, the small and large branch. While it seems challenging to find an analytical 
expression for Q, it is very easy to evaluate it numerically. We have plotted Q, or in other words V2/P 2 , as 
a function of X/P in figure 2. We note that Q decreases monotonically on the large branch, so that, for a 
fixed charge, the A = solution is the most entropic. 




Figure 2: This figure shows V2 as a function of A and P found by numerically solving (3.9). Specifically we plot, 
i>2-P~ 2 = £~ 2 7~ 2 = Q(X/P). The large branch, Ql, is plotted in brown and the small branch, Qs, is plotted in green. 
We should exclude the shaded region, defined by 7 2 £ 2 — 7 2 < 1 , in which, using (3.19), v\ is negative. Since it is 
entirely contained within the shaded region, the small branch is unphysical. The regular black hole, denoted by a blue 
dot, is found on the big branch at X/P — w . At the place where the branches meet, denoted by a red square, v\ — > 00 
(or —00 if we approach from below). 

Having determined V2 (at least in principle), we can find v\ by substituting (3.7) into (3.8) and using 
(3.16), we get, 

vi = f 2 (l-7 _ Ve 2 7 2 -l) _1 (3-18) 
= v 2 (l - [X/P^s/g-^X/P)-!)- 1 , (3.19) 

with (3.12) ensuring reality. 

For v\ to be positive and finite, on sees that from (3.18), we require £ 2 7 2 — 1 > 7 2 - Now, at the branch 
point, the function /(£) = cosh£ — 7^ has a single zero, so we require that /'(£) = when /(£) = 0. In 
other words, in addition to (3.9) the branch point is determined by, 

sinh£ = 7. (3.20) 

Combining (3.9)and (3.20) gives, 

I 2 ? - 1 = 7 2 , (3.21) 

which, using (3.18), implies that as we approach the branch point, v \ — » 00 or in other words the ACLS2 
approaches flat space. From figure 2, we see that the small branch lies entirely in the region 7 2 £ 2 — 7 2 < 1, 
and consequently, v\ is always negative on it, making it unphysical. Discarding the small branch, we have 
used (3.19) to plot V2/V1 as a function of X/P in figure 3. We see that v±/v2 increases monotonically, 
eventually diverging at the branch point. 
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Figure 3: This figure shows v^/vi as a function of X/P found by numerically solving (3.9) and using (3.18). For X/P 
small, we have v\ sa V2- As we approach the branch point, denoted by a red square, v\ diverges and v^/vi — * 0. The 
regular black hole, at X/P — Wq, is denoted by a blue dot. 



Finally, one can actually obtain a large charge, small coupling expansion for Ql. Assuming Ql, has a 
nice Taylor expansion about zero, using (3.17) as a starting point, and by taking successive derivatives of 
(3.13) and (3.15) one can recursively expand Ql(x) about zero. As discussed in appendix A, we find that, 

G L (x) = 1 - x - \x 2 + O (x 3 ) , 
so that we can write a large charge, small coupling expansion for the entropy, 

A 
P 



£ = ttP 2 1 



3 P 2 V P 3 



(3.22) 



(3.23) 



For completeness, we mention that on the small branch, as discussed in appendix A, for x small, we get, 

-2 



g s {x) ^ x[W-i{-\^)] 



(3.24) 



where W_i is the non-principal real branch of the Lambert function. 

As a check we note that on the large branch, taking A — > 0, we recover the usual near horizon extremal 
Einstein-Maxwell Reisner-Nordstrom solution with, 



vi = v 2 



pi 



(3.25) 



4 Entropy function and the regular black hole 

In this section we confirm that the entropy function analysis of the regular black hole reproduces the near 
horizon geometry of the known solution found in [14]. This merely entails considering the results of the 
previous section with the appropriate value of X/P. 

As discussed in section 2, the regular black hole corresponds to the point X/P = Wq, so that (3.9) 
becomes, 

cosh£ = w ~ 1/2 £. (4.1) 
One can analytically check, using the property woe w ° = e" 1 , that (4.1) has a solution, 

or in other words ^{wq) = \{wq + 1). We have plotted the position of the regular solution as a blue dot in 
figures 2 and 3, from which we observe that it is on the large branch. 

Finally, one can check that substituting the solution, (4.2), into (3.16) and (3.18) reproduces (2.16) and 
we are done. 
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5 Conclusion 



In this paper we examine entropy function formalism for regular magnetically charged black hole solution 
in Einstein-Hilbert gravity coupled with a certain non-linear U(l) gauge field. The mass and charge of the 
full solution can be tuned so that it has no curvature singularity at the centre. In the extremal limit this 
corresponds to a particular charge to non-linear coupling ratio with an AdS2 x S 2 near horizon geometry. 
Unsurprisingly we find that the entropy function analysis match with the exact solution found by solving 
the full Einstein equations. This is in contrast with the claim in the recent papers [17, 18]. 

Indeed in the entropy function formalism, the equation of motion, which follow from extremizing the 
entropy function, are exactly the same as equation of motion at horizon found by extremizing the action, 
simply because the entropy function (up to Legendre transformation) is the Lagrangian at the horizon. The 
fact that the entropy is the value of entropy function at its extremum is derived from Wald entropy formula, 
using the near horizon symmetries. Both of these results apparently have nothing to do with the regularity 
of the solution inside the horizon. 

Acknowledgements: We would like to thank Ashoke Sen for helpful comments. The work of K.G. is, in 
part, supported by the EU-RTN network contract MRTN-CT-2004-005104 and INTAS contract 03-51-6346. 
The work of H.Y is supported by the Korea Research Foundation Leading Scientist Grant (R02-2004-000- 
10150-0) and Star Faculty Grant (KRF-2005-084-C00003). 

A Large charge/small coupling expansion of the entropy 

In this appendix we discuss the expansion of Q{x) about zero. 

Taking the derivative of (3.13) with respect to x and solving for T 1 we find, 

T'(x) = coshes)) (A ,x 

{ 1 2^(^sinh(jF(x)) _ l) ' \ A - L > 

while taking of derivative of (3.15) gives, 

g'(x) = -2sech 2 (^(x))tanh(^(x))^'(x) . (A.2) 

Now using (3.13) and (A.l) we can rewrite (A.2) as, 

and using (3.11) and taking the limit x — ► 0, we get, 

g'fQ) = { (small branch) 
\ — 1 (large branch) 

Taking another set of derivative, after some algebra we obtain, 

„ _ ■7 r cosh(2.T) - cosh(^) sinh(^) 
y {X) ~ 2J*(Ftanh(F) - 1)3 ' {A - b) 

and once again taking the x — » limit, we get, 

Q"(0) = { 00 ( sma11 branch) ,^ g . 

\ — | (large branch) 

We found that the second order Taylor expansion of Gl(x), (3.22), agrees well with our numerical plot for 
x small. For example, at x = Wo, we find that they differ by about 2%. 

On the other hand, we see that the small branch does not have a nice Taylor expansion about the origin. 
However, on the small branch, we see from figure 1, that when £ is large, 7 is also large and consequently, 
X/P = 7~ 2 , is small. For £ S> 1, we can approximate (3.9) by, 

ie«« 7 e + 0(e- ? ) , (A.7) 
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which can approximately be solved by, 



(A.8) 



where W is the Lambert function defined in (2.13). For, — e _1 < x < 0, the two real branches of W satisfy 
Wo(x) > —1 and W-i(x) < —1, [20], consequently since we are assuming that £ 3> 1, we should take the 
branch W_i. So, using (3.14,3.15), we obtain, 



which we found agrees well with our numerical results shown in figure 2, for X/P small. 
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